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Single photon sources represent a fundamental building block for optical implementations of quan- 
tum information. In this paper we propose a multiplexed system based on asymmetric configuration 
of multiple heralded single photon sources. The performance and the scalability comparison of the 
scheme to both existing multiple-source architectures and simple faint laser configuration reveals 
an advantage of the proposed scheme in realistic scenarios, and provide insights on the potential of 
these architectures for integrated implementation. 
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I. INTRODUCTION 

Single photon sources (SPSs) represent a key resource 
for optical quantum computing and communication. In- 
deed, optical quantum computers based on integrated 
photonic technology [TH3] can be build by using linear 
optics and SPSs as shown by Knill, Laflamme and Mil- 
burn [5 j. On the other hand, current Quantum Key Dis- 
tribution (QKD) systems [6HS] use weak laser pulses in 
place of single-photon sources and decoy state technique 
[9] [10] to avoid the photon splitting number (PNS) at- 
tack [11] on the pulses containing more than one photon. 
Using true single photon sources would guarantee that 
all pulses contain one and only one photon, increasing 
the key generation rate and avoiding any PNS attack. 

The single photon source usually employed in quantum 
information applications is constituted by a faint laser 
(FL), namely an attenuated and pulsed laser source [23] . 
The number of photons in each pulse can be modeled, for 
a coherent source, by a Poisson random variable K with 
distribution: 

P(K = k) = ^e-», (1) 

where \i is the mean photon number in each pulse, which 
depends on the power of the laser. In order to evaluate 
the output quality of the source two indexes are usually 
employed: the one-photon probability P(K = 1) and the 
signal to noise ratio: 

1-P(K = 0)-P(K = 1) eP-l-n' KJ 

The SNR is the ratio between the one photon probability 
and the probability of having more than one photon in 
the output of the system. This index quantifies a crit- 
ical quantity of the source: multiple photons per pulse, 
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in fact, represent, in optical quantum computing, errors 
whose effects are hard to detect and correct, while, in 
QKD, they provide the possibility of PNS attacks. 

The main limitations of the FL source stems from the 
fact that fi is the only tunable parameter. This induces a 
trade off: the value of \i that maximizes P(K — 1) is given 
by \i = 1 and corresponds to a value of the one photon 
probability of e _1 ~ 0.37. However, for \i = 1 the SNR 
is equal to (e — 2) _1 ~ 1.39, a value that is typically 
unacceptable for applications requiring a single photon 
source. Since the SNR is unbounded for \i approaching 
zero, the mean photon number is usually kept low in or- 
der to avoid multiple photons events, reducing also the 
probability of single photon emissions. To overcome the 
limitations of FL source, in the last years architectures 
with multiple heralded SPS have been proposed. Pho- 
tons are generated by means of spontaneous parametric 
down conversion (SPDC): an intense laser pump imping- 
ing on a nonlinear crystal leads to probabilistic emission 
of entangled pairs of photons into two different spatial 
modes, with their rate depending on the pump intensity. 
It is then possible to "herald" the presence of a photon in 
one of the modes by detecting the correlated twin photon 
in the other one. 

Post-selection schemes in combination with multiple 
heralded SPS have been first proposed by Migdall [12] 
in order to enhance the probability of obtaining a single 
heralded photon, with an implementation that requires 
an m-to-1 global switch. In the same work, the perfor- 
mance of the proposed scheme for finite detection effi- 
ciency were studied. However, as pointed out in [13], an 
efficient implementation of such device is not currently 
available, and it would be hardly scalable. To overcome 
these problems a symmetric scheme, employing a total of 
m — 1 binary polarization-switching photon routers ar- 
ranged in a modular tree structure has been proposed by 
Shapiro and Wong in [13] , where they also evaluate the 
probability of emitting n photons taking into account im- 
perfect detectors and optical switches. An experimental 
implementation of the scheme has been pursued in [14] 
using 4 crystals, along with an essential discussion of its 
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scalability. Recently, the analysis of multimode emission 
in SPDC used for SPS was carried out [T5] . 

However, a thorough analysis of the performance of 
these sources in the presence of finite efficiencies, and 
their comparison with respect to a simple faint-laser 
source, is lacking: nonetheless, it appears to be a crucial 
step in assessing their potential, especially in the light of 
the experimental difficulties reported in [14]. In pursu- 
ing this analysis, we believe that one of the most delicate 
point consists in devising proper performance indexes, 
ensuring a fair comparisons between different methods. 

In this paper, after reviewing the main ideas and the- 
oretical result underlying the existing SPS architectures, 
we propose a new heralding scheme based on an asym- 
metric configuration. Next, we develop a comprehensive 
comparison between the architectures, that investigates 
the comparative performance, the scalability and the lim- 
its of SPS with multiple heralded sources in realistic sce- 
narios. 
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FIG. 1: Schematic of the MHPS scheme. The blue rectan- 
gles labeled with NLC represent the non-linear crystals, the 
detectors are labelled with APD and the orange rectangle, 
labelled with O.S., represents the optical switch. 



II. REVIEW OF HERALDED SOURCES 

A. Multiple sources and 
the advantage of post-selection 

In this section we present the key ideas underlying 
a class of source architectures that outperform the FL 
scheme. The building block of such architectures is rep- 
resented by the so-called Heralded Source (HS). The HS 
exploits the SPDC effect on a non-linear crystal pumped 
with a strong coherent field, which leads (with a cer- 
tain probability) to the simultaneous generation of a 
pair (or more pairs) of photons. One photon of the 
pair is then fed to a photon detector: if the duration 
of the pulse (At p ) is much shorter than the measurement 
time interval (AT), but much greater than the recip- 
rocal of the phase-matching bandwidth i.e., (Agj), i.e. 
AT >> At p » 1/Acj, the statistics of the pairs is still 
Poissonian [15] . 

The HS can be employed in multiple crystal strategies 
that outperform the FL by exploiting the parallel use of 
HS units and post selection strategies: intuitively, the 
advantage of using a scheme exploiting a parallel imple- 
mentation lies in the fact that the intensity of the pump of 
each crystal can be kept low, suppressing thus the multi- 
photons events, while keeping an acceptable production 
rate of single photons. 

Let us assume for now to employ ideal detectors in or- 
der to illustrate in a simple setting the potential advan- 
tages offered by this system. The case of finite efficiency 
will be discussed later in this section. 

We consider standard APD detectors: the information 
that these detectors provide is limited, since they are only 
able to discriminate between the case of no incident pho- 
tons and the case in which photons are detected, without 
resolving their number. When a detector is hit by the 
photons, it returns an electric signal (trigger), which in- 



dicates the presence of at least a photon in the channel. 
A trigger event may be described by a Bernoulli random 
variable r with values or 1. 

Let us consider a Multiple Crystals Heralded Source 
Architecture with Post-Selection (MHPS) built as fol- 
lows: an array of m HS units, labelled with index 
i = l,...,m, and each one simultaneously fed with a 
laser pulse with intensity such that the mean number of 
produced pairs is ft. For each HS unit, one mode, cor- 
responding to one of the photons in the SPDC pair, is 
used as trigger and the other is injected into an opti- 
cal switch. Let and Ni be the trigger variable and 
the random variable (r.v.) accounting for the number of 
photons generated at the z-th source. We introduce a r.v. 
X, called post-selection function, that is a deterministic 
function of ti, . . . , r m such that: 

X(ri,...,r m ) e {0,1,..., m}. (3) 

The outcome of x represents the channel that we se- 
lect to be routed to our global output: this channel will 
correspond to the output of the z-th source when x £ 
{l,...,m}, otherwise, when x = 0, no channel will be 
selected. We then denote N = N x as the r.v. accounting 
for the output statistics. 

In other words, we assume to have an optical switch 
that is able to select one channel among the others ac- 
cording to the selection rule specified in (J3|. As first 
proposed in [12], we may assume: 

X(n, . . . , r m ) = z with r* = 1 and 75 <» = . (4) 

Furthermore, we assume that if all the detectors do no 
fire: 

X(r 1 = 0,...,r m =0) = 0. (5) 

This choice of the the x function correspond to the fol- 
lowing strategy: the output signal is taken from the first 
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source (starting from i = 1) that triggers (thus indicating 
the presence of at least a photon in the channel). It is 
worth noting that the precise structure of x is n °t impor- 
tant in this ideal situation. In fact, any selection function 
satisfying x(ri,...,r m ) = i if at least n = 1 could be 
used in place of Q without altering the performances. 
We will see in the next section how different selection 
rules affects crucially the performances once non-ideal 
situation are kept into account. 

In order to derive the output statistics for the MHPS, 
let us first compute the probability that \ 7^ : 



P(X + 0) = 1 - [] P(t 4 = 0) = 1 - e 



(6) 



as well as the probability of having n photons in the 
global output, given that x — i'- 



P(N = n\x = i) 



Vjil-Sn), ifi^O 

if i = 0. 



(?) 



In the previous expression S n is the Kronecker Delta 
(So — 1 and 5 n ^o = 0). We thus have that the prob- 
ability of having n photon in the global output of the 
MHPS is given by: 



?(N = n) =?(N = n\ X ^ 0)P(* ^ 0) + 
+ P(7V = n| X = 0)P(x = 0) 



(8) 



1 



In particular, for the N = 1 case: 
P(N = 1) =Jle- J1 — 



^{i-5 n ) + 5 n e- m ». 



1 - e~v 



(9) 



Let us now compare the performances of the MHPS 
with the performances of FL. Using ([2| we find that the 
expression of the SNR for the MHPS is given by: 



SNR 



lie 



fie 



(10) 



Notice that, when jl = /i, the signal to noise ratio is equal 
to the SNR of the Faint Laser, while the single photon 
probability P(N = 1) is always larger. In Table [i] we 
listed the one photon probability for different values of 
the SNR and of the number of crystals. The benefits of 
the MHPS with respect to the FL are apparent: fixed 
any SNR level, it is possible to obtain a higher value of 
the one photon probability. This is because of the post 
selection procedure, that can turn (with a certain proba- 
bility) events in which more than one detectors trigger at 
the same time into an event that corresponds to a one- 
photon output by blocking the output of all the HS units 
but one. 



TABLE I: One photon probability for the FL and the MHPS 
scheme, with m = 4, 8 and 16. 
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FIG. 2: Schematic of the binary switch. 



Before closing this section we consider the ideal per- 
formances of the MHPS performances in the limit of an 
infinite number of crystals. In this case, in order to avoid 
the need of infinite pump power, we must rescale the in- 
tensity as ji — » ji/m. This is equivalent of using a parallel 
strategy in which a single laser beam of intensities \i is 
split into m beams with intensity jl = fi/m shining m 
different crystals. In the infinite number of crystal limit 
we get: 



lim P(N = 1) = 1 - e" 



(11) 



This expression can be rewritten as 1 — (e~^) m , repre- 
senting the complement of the probability of not having 
triggers among all the detectors, i.e, the probability that 
at least one among them triggers. This means that, since 
we are using an infinite number of crystals, the intensity 
per single crystal goes to zero, and so does the probability 
of generating multiple photons in a single crystal: each 
crystal is going to produce either one or no pairs. This 
also implies that the MHPS achieves an infinite value for 
the SNR as can be easily seen from (10). 



Finite efficiency and symmetric modular 
architecture 



Any discussion regarding the physical implementations 
would be vain without taking into account the realistic 
efficiencies in detection and routing of the produced pho- 
tons. An actual detector, in fact, is subject to losses 
whose effects are usually described by introducing a pa- 
rameter < 7] < 1, called detection efficiency, that repre- 
sents the probability of detection of an incident photon. 
This parameter takes into account also the collection effi- 
ciency, as the phase matching relations yield uncertainty 
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TABLE II: Definition of the selection function \ 



in the direction of the emitted photons. 

Physical implementation of post-selection rules are 
subject to losses as well. The efficiency in transmission is 
modeled with a parameter < 7 < 1, called transmissiv- 
ity, that represents the probability of transmission for a 
photon through the router. It is important to note that 
this limited efficiency is referred only to the transmission 
of the photons: for what it concerns the transmission of 
the electrical signal we are always going to assume that, 
once a trigger happen, it is transmitted until the end of 
the transmission chain without errors. 

We remark that, in order to consider the role of a 
finite transmissivity, it is key to specify the particular 
routing/switching architecture that is being considered, 
since the potential gain with respect to a FL (without 
the routing inefficiencies) will in general depend on it. 
This is not the case in the ideal scheme described in the 
previous section: expressions 



and (10) do not depend 
on the particular x chosen as long as it satisfies Q and 

©. 

In what follows we are going to briefly review a mod- 
ular architecture, proposed in [13 , that employs binary 
photon routers (2-to-l) whose working mechanism is de- 
picted in fig. |5J tr^l £ {0, 1} are Bernoulli r.v. account- 
ing for the trigger signal incoming from the previous level. 
We introduce the r.v. x, called selection function, that 
is a deterministic function of tl and tr as defined in ta- 
ble |n[ Therefore, x accounts for which channel is driven 
through the router. 

The m-HS units are arranged as shown in fig. [3j where 
r-j, with i = 1, . . . , m, are r.v. accounting for the trigger 
of the i-th channel. The outputs of the first stage are fed 
into the second stage's routers and so on, until the end of 
the transmission chain. This architecture can be clearly 
realized only for m = 2 fc , with k G N. 

It is worth noting that any successfully transmitted 
photons have to pass trough k = log 2 m routers. We will 
call this scheme Symmetric Multiple-crystals Heralded- 
Source with Post-selection (SMHPS). The overall effect 
of the selection function implemented in the routers is 
that, if more than one detectors fire, the channel routed 
to the end of the chain is the one coming from the crys- 
tal corresponding to the lowest i. More precisely, this 
architecture is efficiently implementing the post-selection 
function x given in Q. 

The corresponding probability of having n photons in 
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FIG. 3: Schematic of the SMHPS scheme proposed in [T3] . 
The blue rectangles labeled with NLC represent the non- 
linear crystals, the detectors are labelled with 77 and the or- 
ange squares, labelled with 7, are the photon routers. 



the final output can now be computed to be: 

P(N = n) = e~ r]2k v5 n + 
(/i 7 /c ) n e~^ k 1 - (1 - ^V^ 1 ^)^ 



1 - 



-(1 



(12) 

) 



The first term in the previous sum account for the prob- 
ability of having no photons in the output because no 
detector fired. We will postpone the performances com- 
parison of SMHPS with the FL in section [TTTj 



III. ASYMMETRIC ARCHITECTURE 

A. Analysis of the 2-crystal architecture: 
Symmetry is not optimal 

In the all the previous works, it was assumed that all 
the crystals were driven with same intensity. We will 
now prove that this choice represents a suboptimal case 
for the one photon probability. 

In order to do so, let us derive the probability of having 
one photon in the output of a MHPS composed by two 
crystals each one be fed with different intensities such 
that the mean number of emitted pairs is fij . Let us call 
N the r.v. accounting for the statistics of the final output 
while N\ (N2) and t\ (t<i) accounts for the statistics and 
trigger of the first (second) channel respectively. Let x 
be a post-selection function as specified in the previous 
section. 

The probability that the first and the second detector 
give rise to a trigger is given respectively by: 



P( X =1) = (1 
P(X = 2) = (1 



-Mi), 
-M) e - 



(13) 
(14) 
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FIG. 4: Contours of the one-photon probability (15) for the 
two crystals MHPS. The blue line represent the \i\ — \±2 case. 



FIG. 5: Schematics of the AMHPS. Notice that each crystal 
is fed with different intensities to compensate the different ab- 
sorption rate of different channels. The blue rectangles labeled 
with NLC represent the non-linear crystals, the detectors are 
labelled with 77 and the orange squares, labelled with 7, are 
the photon routers. 



while the probability of having one photon in the final 
output is given by: 



P(N = 1) = me 



mi 



/i 2 e 



-O1+M2) 



(15) 



The maximum of P(N = 1) is achieved when 

(^ 2 ) = (l-e-\l). (16) 

This can be seen in fig. [4] where the contours of P(N = 1) 
are plotted: the blue line represent the cases \±2 = fii- 

Therefore, we conclude that, in general, feeding all the 
crystals with the same laser intensities, leads to subopti- 
mal performances. It is worth noticing that the asymme- 
try comes from the selection function x of Table [3j where 
the first source is initially checked and, only if this source 
doesn't trigger, the selection function consider the second 
source. As we demonstrated, this asymmetry can be ex- 
ploited to improve the probability of having one photon 
at the output channel. 

In the ideal case (perfect detection and transmission), 
it is straightforward to compute the intensities that lead 
to the one photon probability optimal value generalizing 
the analysis above. Anyway, such optimization task be- 
come quite non-trivial if finite efficiency are taken into 
account, especially for a large number of crystals. 



B. Proposed asymmetric architecture 

Here we would like to propose an alternative (still sub- 
optimal) scheme which can be realized with an arbitrary 



number of crystals and achieves better performances than 
the SMHPS in many situation of experimental interest. 

Let us again suppose to have an array of m-HS system. 
Let Ti and Ni be the r.v. accounting for the trigger and 
the number of photons in the input of the i-th channel 
respectively. Finally \±i is the mean number of generated 
pair of the i-th crystal. This new scheme also employes 
binary switches but the multiplexing is performed in a 
different way, as displayed in fig. [5] the m-th source is 
multiplexed with the (m-l)-th, the output is then mul- 
tiplexed with the (m-2)-th source, and so on. In this 
new configuration each successfully transmitted photon 
passes through a number of photon routers ki that de- 
pends on which crystal it has been created from. More 
precisely, it is given by: 



ki 



if i ^ m - 
1, if i = m. 



1 



(17) 



Since each channel is subjected to a different attenuation 
(see (17)), we compensate the different absorption rates 
by choosing the input intensities as follows: 



Mi 



r 



k i 



l,...,m. 



(18) 



This choice is still suboptimal but is sufficient to out- 
perform the SMHPS in many situation of experimental 
interest. 

This multiplexing architecture is therefore asymmetric 
and will be denoted with AMHPS, to highlight the dif- 
ferences with the SMHPS scheme. It is worth to note 
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that both the symmetric than the asymmetric architec- 
ture require the same number of detectors and routers. 
With each binary switch configured as before, the chan- 
nel that we select to be routed to our global output, is 
here determined by: 



AMHPS 



x 



1 and Tj 



for j = 1, . . . ,i - 1, 



if Tj = for j = 1, . . . , m. 



(19) 



It is worth noting that, the global effect of this rule, 
within this architecture, is that if two or more channels 
are heralded, the one that needs to pass through less 
routers is selected and routed to the end of the chain, 
thereby decreasing the probability of absorption. 

As shown in appendix [Aj the probability of emitting 
n-photons for the AMHPS is given by: 



P(N = n) = 5 n e-^ 

m 



(2-7)7 ± 



(20) 



1 - (1 - r]) n e^e 



It worth noting that, by using the compensation pro- 
posed in (18), the dependence of the one photon prob- 



ability on the intensities is reduced to a single variable, 
i.e. ft. Moreover, when 7—^1, the value of P(N = n) for 
the asymmetric scheme coincide with the values (12) of 
the symmetric scheme for any value of 77. 

In the next section we are going to confront the FL, 
the SMHPS and the AMHPS with the compensation we 
have proposed. 



IV. PERFORMANCE COMPARISON 




Log 2 m 



FIG. 6: One-photon probability for the AMHPS with 6 = 10 
and m £ {4, . . . , 256} and various pairs of (77, 7). 
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One-photon probability for the SMHPS with = 10 
G {4, . . . , 256} and various pairs of (77, 7). 



A. Maximum probability with guaranteed SNR 

This section is devoted to the comparison of the FL, 
the SMHPS and the AMHPS by means of numerical con- 
siderations. As we have stated previously, in many appli- 
cation it is crucial to be able to rely on a threshold value 
for the SNR. With this in mind, we propose the following 
method for the comparison: we fix a threshold value for 
the acceptable SNR, O, and for each architecture, given 
pair of efficiencies (77, 7) and for a given number of crys- 
tals, we compute the optimal pump intensity for which: 

fifa 7, 6) = arg max P(7V = 1, /i, 77, 7), (21) 

i.e., we compute the intensity that maximize the one pho- 
ton probability provided that the SNR has a greater or 
equal value than 9. We notice that, since lim^o SNR = 
+00 and lim /7 ^ +oc SNR = for both the symmetric 
than the asymmetric scheme, by choosing the appropri- 
ate value of /i, any value of the SNR can be achieved. 

We finally compare the one photon probability for dif- 
ferent methods constrained to the optimal intensities: 

P(JV = l) = P(JV=l,/i(»7,7,©),'7,7)- (22) 



Notice that, once the number of crystals and the SNR 
threshold are fixed, the above quantity depends only on 
(77,7). From now on, P(N = 1) always indicates this 
optimized probability with the SNR constraint. 



B. Scalability of the schemes for finite efficiencies 

We first discuss the scalability features of the strategies 
with respect to the total number of crystals. In fig. [6} 
we plotted the values of P A (N = 1) and F 5 (7V = 1) 
versus the number of crystals, m, ranging from 4 to 256 
for different pairs of (77, 7) and for threshold SNR given 
by 6 = 10. 

For what concern the AMHPS it is apparent that, for 
all the considered pairs of (77,7), the one photon proba- 
bility increases until it reaches an asymptotic value (the 
dependance of this value on 7 and 77 is non-trivial). This 
fact implies that, once the detection efficiency and the 
transmissivity are fixed, there is a threshold value on 
the number of crystals above which there is no further 
improvement in the performances of the scheme. It is 
worth noting that for experimental realistic parameters, 
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i- e - 7 ^ 0-5, the asymptotic performances are practically 
already reached for m = 8. 

For what it concerns the SMHPS, after an initial tran- 
sient, the one photon probability starts to decrease (a 
part for the ideal case of 7 = 1). We have analytically 
shown that, excluding the 7 = 1 case, in the limit of 
infinite number of crystals (m —> +00) the one photon 
probability of the symmetric scheme approach the one 
photon probability of the faint laser when the pump pa- 
rameters are chosen in order to have asymptotically the 
same SNR (as can be sees in fig. [7]). The result is demon- 
strated in appendix |C] 

Let us try to give a motivation for this counter in- 
tuitive behavior: has we have mention above, in the 
SMHPS each successfully transmitted photon have to 
pass through log 2 m routers, increase the number of crys- 
tal means also to increase the absorption rate the photons 
are subjected to. Thus for this geometry architecture 
the benefits deriving from the increase of the number of 
crystals do not compensate the increase in the absorp- 
tion rate. On the contrary, as we have mention above, in 
the AMHPS are most likely to be selected those channels 
whose photons have to pass through less routers in order 
to reach the global output leading, on average, to a lower 
absorption rate. 

Summarizing, there are some benefits for the SMHPS 
in increasing the number of crystals but only up to a cer- 
tain number, depending on the detection efficiency and 
transmissivity. Anyway increasing further the number of 
crystals will lead to poorer and poorer performances. 

The AMHPS offers significant benefits in increase the 
number of crystal until a certain number depending 
on the detection efficiency and transmissivity, once the 
threshold is reached increasing further the number of 
crystal will left the performances unchanged. Finally, 
we remark that since the two methods employ the same 
post-selection rules, the gap in the performances arises 
form the different architecture geometries that is respon- 
sible for the different distribution of the routers. 

Before closing this section, let us briefly review and 
discuss the scalability analysis proposed in [12]. In that 
paper, in order to evaluate the advantages of a scheme 
with respect to the FL, the gain, namely the ratio G 
between the one-photon probability of having one pho- 
ton in the output of the SMHPS and the probability of 
producing one photon with the FL, is considered. It is 
worth remarking that in [12] the one-photon probability 
for the proposed scheme is computed neglecting both de- 
tection and transmission inefficiencies: in this ideal case, 
it turns out that both the scheme and the FL have the 
same SNR provided that the intensity with which the 
HS units are fed is equal to the intensity of the FL. In 
order to analyze the scalability taking into account the 
absorption due to the routing chain, they propose to con- 
sider the asymptotic behavior of the product j k G, thus 
comparing the benefits of a growing multi-crystal archi- 
tecture to the increase in the absorption rate due to the 
longer routing chain. As a result of this analysis, we 



have that the advantage of the SMHPS is maintained 
(i.e. lim^ oo7 fe G>l) if 7 > 1/2. 

In order to perform a similar analysis we should com- 
pare the asymptotic behavior, in the limit of an infinite 
number of crystals, of the rate between the one-photon 
probability for the SMHPS, (p), with the probability of 



producing one photon with the FL, for the same SNR. 
As shown in appendix [Cj in the k — » 00 limit the SNR 
of the faint laser with intensity /i is equal to the SNR 
of the SMHPS with fx = (the number of crystal is 

m = 2 k ). Moreover, as explained previously, in order to 
avoid infinite power we need to rescale the pump power 
(for both the FL than the SMHPS) with the number of 
crystal, /i — » fi/2 k . The gain we obtain is given by: 



G 



(1 



')[1 



(-1+7*) 



(27) fc (rq - 1)] 



1 



" (2 7 )fc 



(23) 



In the infinite k limit, the gain G is greater than 1 for 
7 > \ (and is actually divergent for 7 > 1/2) as found in 
[T2] . However, in the rescaled pump power case, the one- 
photon probability of the SMHPS tends asymptotically 
to zero and the gain does not seems the proper index to 
evaluate the absolute performance of a scheme. In fact, 
despite exhibiting an advantage with respect to the FL, 
P(7V = 1) is asymptotically zero in both cases (as already 
said, without the rescaling \i — > /i/2 fc , the single photon 
probability of the SMHPS and of the FL coincide in the 
large k limit, and the gain is always 1). This means that, 
with finite transmissivity, not only increasing the num- 
ber of crystal beyond a certain value does not bring any 
advantage, but it is actually detrimental to the SMHPS 
scheme performance. On the other hand, the AMHPS 
scheme is "robust" with respect to implementations with 
large numbers of crystals, as it is clearly shown in fig. [6] 



C. Comparison between AMHPS and SMHPS 

In figures [8] and [9] we displayed the contours plot of 
P(N = 1) for the AMHPS and the SMHPS with SNR 
threshold given by 9 = 10 and m = 4 and 8 respec- 
tively. The white zones correspond to the case in which 
P(N = 1) is below the value of the one-photon probabil- 
ity of the FL corresponding to an SNR equal to 10 (for 
the faint laser, when SNR=10 the one photon probabil- 
ity is given by 0.155). We can see that both the AMHPS 
and the SMHPS outperform the FL in a vast region of 
the plane (77, 7), but such region is wider for the AMHPS. 
As expected, the best performances are reached for high 
values of both the detection efficiency and the transmis- 
sivity: furthermore in this limit the two methods are 
comparable since for n — >• 1 and 7 — >• 1 they both tend to 
the ideal situation of the MHPS. 

In the light of the numerical results, we can conclude 
that systems that exploit post-selection can indeed of- 
fer benefits with respect to the FL, for a wide set of 
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FIG. 8: Contours of the one-photon probability for the SMHPS (left) and the AMHPS (right) with guaranteed SNR, 6 = 10, 
and m — 4. In the white area P(N = 1) is below the value of the one-photon probability (given by 0.155) of the FL with SNR 
equal to 10. 
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FIG. 9: Contours of the one-photon probability for the SMHPS (left) and the AMHPS (right) with guaranteed SNR, 9 = 10, 
and m = 8. In the white area P(N = 1) is below the value of the one-photon probability of the FL with SNR equal to 10. 



values of experimentally achievable detection efficiencies 
and transmissivities. 



Let us now compare the AMHPS and the SMHPS. 
Let as call P A (N = 1) and P S (N = 1) the one-photon 
probability as computed in (22) for the AMHPS and the 



SMHPS respectively and let us define: 

^ A {N = 1)-P S (N=1) 



A := 



Ps(N = l) 



(24) 



the percentage differences between the two single photon 
probabilities. 



In fig. 10 and 11 we shown the contour of A for an 
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FIG. 11: Contours of the percentage differences (24) with guaranteed SNR, G = 10, and m = 16 (left) or m = 32 (right). A 
positive number correspond to an advantage of the AMHPS to the SMHPS. In the white area Pa(N = 1) and Ps(iV = 1) are 
both below the value of the one-photon probability of the FL with SNR equal to 10. 
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SNR equal to 10 and m = 4,8 i 16 and 32. The white 
zone are those in which neither Pa(N = 1) nor Ps(N = 
1) are grater than the FL one-photon probability, that 
for and SNR of 10 is equal to 0.156. 

In the case of m = 4 and m = 8 the AMHPS always 
outperforms the SMHPS. In the case of m = 16 and 
m = 32 the AMHPS outperforms the SMHPS in a vast 
portion of the plane (77, 7) even if a zone in which the 
SMHPS slightly outperforms (less than 2% for m = 8 
crystals) the AMHPS is present for high value of 7. 
Anyway, the advantage of the AMHPS in more realis- 
tic situations is apparent, especially in the area where 
7 ^ 0.5 and the detection efficiency is higher than 0.5, 
the AMHPS significantly outperforms the SMHPS. Any- 
way, above m = 32 the AMHPS outperforms the SMHPS 
in every portion of the plane (77, 7) again. 



V. EXPERIMENTAL FEASIBILITY AND 
CONCLUSIONS 

Let us now discuss a possible experimental realization 
of the proposed SPS configuration. Nowadays, integrated 
devices represent the best resource to achieve high effi- 
ciency of the SPDC process and ensure good coupling 
into single mode fibers (for a review on integrated source 
see [H]). It is possible to use non-degenerate collinear 
phase matching and a dichroic beam splitter (or alterna- 
tively using a counterpropagating mode source [17]) in 
order to separate the two beams. For instance, it was 
recently reported the possibility of heralding single tele- 
com photons at 4.4 MHz rate with 45% heralding effi- 
ciency [18]. Moreover, to efficiently detect the triggered 
photon, high efficient transition-edge sensors (TES) can 
be used: an heralding efficiency of 77 ~ 62% has been 
recently reported by using TES [19] (see also [20 for a 
review on single photon detectors). Regarding the op- 
tical switch, a 2x2 silicon electro-optic switches with a 
broad bandwidth (60 nm), an ultrafast speed (6 ns) and 
a transmission of 7 ~ 50% has been reported [21]; other 
modulators otherwise allow lower losses at the cost of re- 
duced working spectrum [22 . As shown in Fig. [8] with 
these values of 77 and 7 the heralded schemes outperform 
the faint laser even with low number of crystal (m = 4), 
with the asymmetric scheme more performant than the 
symmetric one. 

In conclusion we have proposed an asymmetric archi- 
tecture for the multiplexed heralded single photon source 
and we have compared it with the symmetric version pro- 
posed in [13 and with the faint laser source. We have 
proven that the asymmetric architecture outperform the 
FL and the symmetric scheme in a vast region of the pa- 
rameter space (77,7). We have also demonstrated that, 
in the large number of crystal limit and by considering 
a fixed SNR, the symmetric configuration is asymptot- 
ically equivalent to the faint laser for any 7 / 1, while 
for the asymmetric scheme the one photon probability in- 
creases until it reaches an asymptotic value dependent on 



7 and 77. Values close to the asymptotic ones, at least for 
experimentally available efficiencies, are reached already 
around 8 crystals. This implies that, even if the expected 
energetic consumption for the AMHPS scheme is higher 
due to the implemented pre-compensation for the losses 
in the routing chain, the necessary overhead will be very 
limited. We believe that our results will be relevant to 
any future realization of heralded single photon source 
based on multiplexed architecture. 
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Appendix A: Statistics of the SMHPS 

For the sake of completeness, let us provide a deriva- 
tion for the SMHPS statistics presented in [13]. Given 
the post selection function x of Q, the probability that 
X = i is given by: 

vl Jp(ri = i) nti p(-n = 0), if^o, 
ln?=iPfa = o) if* = o. 



'(l_ e -W) e -(i-l)W ; ifi^Q, 



if i = 0. 



(Al) 



The probability of having n photons in the i-th. channel, 
given that x = i'- 

{^r-e _/i (l — (1 — r)) j ) 
' iH ^ (A2) 
Sj, if z = 0. 

The probability of having n photons in the final output 
given that x = * an d Ni = j ^ n is given by: 

P(N = n\ X = i,Ni=j)= (A3) 

= f (D(VT(l-7 fc ) j ~ n , ifi^Oandn^ j 

\Sn, if 2 = 0. 

Finally the probability of having n photons in the final 
output is given by: 

2 k 00 

P(N = n) = E P ^ = n\ X = i,N z = j) x (A4) 

2=0 j—n 

x P(JVi=.7| X = *)P(x = *) 
=5 n e- ,l2k » + [1 _ (1 _ ^"e-^i-T^A" x 



2" 

i=i 



(i-DvH (At )"e~ 
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For which we have finally: 

P(N = n) ^e-^^Sn + (1 - e-^)x 

(fl~i k ) n e~^ k 1 - (1 - v )n e -ri(l-l k )» 



(A5) 



1 - e-w 



Appendix B: Statistics of the AMHPS 



By computing the products we obtain (20). 



Here we explicitely derive eq. (20). Let N denote the 



r.v. accounting for the output statistics. The probability 
that a trigger coming from the i-th detector is routed to 
the global output is given by: 



>(^ = i)n^iP(^ = o), if <^o, 

ntiHri = 0) if t = 0. 

'(1 - e"^*) nj=J if i ^ 0, 

Il^i e"^ 1 if z = 0. 



(Bl) 



The probability of having n photons in the i-th chan- 
nel, given that \ — ^ 



Appendix C: Infinite crystals limit of the SMHPS 



We here analytically show that the one photon proba- 
bility of the symmetric scheme approaches the one pho- 
ton probability of the faint laser when the pump inten- 
sities are chosen in order to have the same SNR. Let's 
consider the SHMPS with jl = where \i is the faint 

laser intensity and 2 k are the number of crystals. This 
is the natural choice in order to compensate the effect 
of the k routers that each emitted photon has to travel 
through. In this case the SNR of the SHMPS becomes 



P(Ni=j\x = i) 



i(l-(l-rjy) 



(l-e- 







Sj, 



if i ^ o (B2) 
if i = 0. 



The probability of having n photons in the final output 
given that x = * an d Ni = j ^ n is given by: 

P(N = n\ X = hN l =j)= (B3) 

_ f (i)h ki ) n (l - j ki ) j - n , if z ^ and n < j 
\Sn, if i = 0. 

It is worth noting that, the requirement that n ^ j is 
due to the fact that the routing chain can only absorbs 
photons. Finally the probability of having n photons in 
the final output is given by: 



e t e 



k p rjfi 



(1-7?) 



[e^ + e^(ry/i-/i- 1)] 



(CI) 



In the large k limit, when 7 ^ 1, the previous expression 

For 



is equal to the SNR of the faint laser, namely 
this choice of the pump parameter, namely jl = the 
SNRs are equal in the large k limit and the one photon 
probability of the SHMPS becomes 



P(N = n) = J2J2 P ^ = n l* = i > N i= J) ; 

i=0 j=n 

x P(Ni=j\ x = i)P(x = i) 

/ m \ m / i—1 



=s n m< 



\i=i 



eh 

i=l \Z=l 



(Ai7 fe< ) n e" /Ii7 



(B4) 



P(iV = 1) 



He-^[l - e" ( ^ ) - e ^ e^(l - 77)] 



_ ZZii 

1 — e ^ k 



fie 



for k — » 00 and 7^1, 



(C2) 



asymptotic, in the large number of crystal limit, to the 
one photon probability of the faint laser. 
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